Exceptional Point Dynamics in Photonic Honeycomb Lattices with VT Symmetry 

Hamidreza Ramezani and Tsampikos Kottos 
Department of Physics, Wesleyan University, Middletown, Connecticut 06459, USA and 
MPI for Dynamics and Self-Organization - Bunsenstrafie 10, D-37073 Gottingen, Germany 

Vassilios Kovanis 

Air Force Research Laboratory, Sensors Directorate, Wright Patterson AFB, OH 45433 USA 

Demetrios N. Christodoulides 
Gollege of Optics & Photomcs-CREOL, University of Gentral Florida, Orlando, Florida 32816, USA 

(Dated: December 21, 2011) 

We theoretically investigate the flow of electromagnetic waves in complex honeycomb photonic 
lattices with local VT symmetries. Such VT structure is introduced via a judicious arrangement 
of gain or loss across the honeycomb lattice, characterized by a gain/loss parameter 7. We found a 
new class of conical diffraction phenomena where the formed cone is brighter and travels along the 
lattice with a transverse speed proportional to y/^. 



PACS numbers: 42.25.Bs, 11.30.Er, 42.70.Qs 

I. INTRODUCTION 

Conical refraction phenomena i.e the spreading into a 
hollow cone of an unpolarized light beam entering a bi- 
axial crystal along its optic axis, is fundamental in classi- 
cal optics and in mathematical physics [iHl]- Originally 
predicted by Hamilton in 1837 [3] and experimentally 
observed by Lloyd (4j , these phenomena have been inten- 
sively studied in recent years, by a large community of 
theorists and experimentalists [THS]. The physical origin 
of the phenomenon is associated with the existence of the 
legendary diabolical points, that emerge along the axis of 
intersection of the two shells associated with the wave 
surface. Around a diabolical point the energy dispersion 
relation is linear while the direction of the group velocity 
is not uniquely defined. Recently conical diffraction was 
observed in 2D photonic honeycomb lattices [5j which 
shares key common features, including the existence of 
diabolical points, with the band structure of graphene 
in condensed matter physics literature. In graphene, the 
electrons around the diabolical points of the band struc- 
ture behave as massless relativistic fermions, thus result- 
ing in extremely high electron mobility. Both photonic 
and electronic graphene structures allow us to test ex- 
perimentally various legendary predictions of relativistic 
quantum mechanics such as the Klein paradox and 
the dynamics of optical tachyons [TU] . 

While diabolical points are spectral singularities as- 
sociated with Hermitian systems, for pseudo-Hermitian 
Hamiltonians, like those used for the theoretical descrip- 
tion of non- Hermitian optics, a topologically different sin- 
gularity may appear: an exceptional point (EP), where 
not only the eigenvalues but also the associated eigen- 
states coalesce. Pseudo-Hermitian optics is a rapidly de- 
veloped field which aims via a judicious design that in- 
volves the combination of delicately balanced amplifica- 
tion and absorption regions together with the modulation 



of the index of refraction, to achieve new classes of syn- 
thetic meta-materials that can give rise to altogether new 
physical behavior and novel functionality [11] [T2] . The 
idea can be carried out via index-guided geometries with 
special antilinear symmetries. Adopting a Schrodinger 
language that is applicable in the paraxial approxima- 
tion, the effective Hamiltonian that governs the optical 
beam evolution is non -Hermitian and commutes with the 
combined parity (V) and time (T) operator [I3HI5] . In 
optics, VT symmetry demands that the complex refrac- 
tive index obeys the condition n(r) = n*(— r). It can be 
shown that for such structures, a real propagation con- 
stant (eigenenergies in the Hamiltonian language) exists 
for some range (the so-called exact phase) of the gain/loss 
coefficient. For larger values of this coefficient the sys- 
tem undergoes a spontaneous symmetry breaking, corre- 
sponding to a transition from real to complex spectra 
(the so-called broken phase). The phase transition point, 
shows all the characteristics of an exceptional point (EP) 
singularity. 

VT symmetries are not only novel mathematical cu- 
riosities. In a series of recent experimental papers VT 
dynamics have been investigated and key predictions 
have been confirmed and demonstrated |16H19j . Symme- 
try breaking has been experimentally observed in non- 
Hermitian structures |16H18| while power law growth 
- characteristic of phase transitions - of the total en- 
ergy has been demonstrated close to the exceptional 
points in Ref. |18j . In a silicon platform claims have 
been made that non-reciprocal light propagation in a 
silicon photonic circuit has been recorded [TS]. VT- 
synthetic materials can exhibit several intriguing fea- 
tures. These include among others, power oscillations 
and non-reciprocity of light propagation [TTl [1^1 [201 HI] , 
non- reciprocal Bloch oscillations and unidirectional 
invisibility [23^ . More specifically, a recent paper has pro- 
posed photonic honeycomb lattices with VT symmetry 
[TU] . Interestingly, that work has shown that introducing 



alternating gain/loss to a honeycomb system prohibits 
the T^T-symmetry, but adding appropriate strain (di- 
rection and strength) restores the symmetry, giving rise 
to T^T-symmetric photonic lattices [TU]. Moreover, the 
mentioned work has found that in such systems, much 
higher group velocities can be achieved (compared with 
non PT-symmetry breaking systems), corresponding to 
a tachyonic dispersion relation (lOj . In the nonlinear do- 
main, such pseudo-Hermitian non-reciprocal effects can 
be used to realize a new generation of on-chip isolators 
and circulators [21]. Other results within the framework 
of T'T-optics include the realization of coherent perfect 
laser absorber [55] , spatial optical switches [21] , and non- 
linear switching structures [22] ■ Despite the wealth of 
results on transport properties of T-'T-symmetric ID op- 
tical structures, the properties of high dimensional VT 
optical lattices, (with the exception of few recent studies 
[HI EH]), has remained so far essentially unexplored. 

Recently, it was pointed [29] that PT-symmetric 
Hamiltonians are a special case of pseudo-Hermitian 
Hamiltonians i.e. Hamiltonians that have an antilinear 
symmetry |30H32) . Such Hamiltonians commute with an 
antilinear operator ST, where 5 is a generic linear opera- 
tor. The corresponding Hamiltonian H is termed general- 
ized 7^7~-symmetric. In a similar manner as in the case of 
T'T-symmetry, one finds that if the eigenstates of H are 
also eigenstates of the iST-operator then all the eigenval- 
ues of H are strictly real and the iST-symmetry is said to 
be exact. Otherwise the symmetry is said to be broken. 
An example case of generalized PT-symmetric optical 
structure is the one reported in [33] (see also [20]). The 
unifying feature of these systems is that they are built 
of a particular kind of "building blocks" which below are 
referred to as dimers. Each dimer in itself does have VT- 
symmetry and it can be represented as a pair of sites with 
assigned energies e„ , e* . The system is composed of such 
dimers, coupled in some way. For an arbitrary choice 
of the site energies e„ and coupling between the dimers, 
the system as a whole does not possess PT-symmetry 
(indeed, such global T^T-symmetry would require precise 
relation between various e„ and coupling symmetry be- 
tween the dimers). On the other hand, since each dimer 
is T'T-symmetric, with respect to its own center, there 
is some kind of "local" T^T-symmetry (which we shall 
define as T^^T-symmetry) . The main message of these 
papers is that, T^^T-symmetry ensures a robust region 
of parameters in which the system has an entirely real 
spectrum. Below, we will use the term 'PT'-symmetry 
in a loose manner and we will include also systems with 
generalized T'T-symmetry. 

In this paper we investigate beam propagation in 
non-Hermitian two dimensional photonic honeycomb lat- 
tices with T'T-symmetry and probe for the possibility 
of abnormal diffraction. We find a new type of conical 
diffraction that is associated with the spontaneously VT- 
symmetry breaking phase transition point. Despite the 
fact that at the EP the Hilbert space collapses the emerg- 
ing cone is brighter and propagates with a transverse ve- 





FIG. 1: Honeycomb photonic lattice structure with intra- 
dimer coupling t and inter-dimer coupling ta = 1. Sub- lattice 
(lossy waveguide) a„,„i is shown by green circles while sub- 
lattice (gain waveguide) fe„,m is shown by the red circles. Each 
dimer is distinguished by index n and m. The field is coupled 
evanescently between the waveguides. 



locity that is controlled by the gain and loss parameter 
7- 

The organization of the paper is as follows: In Section 
In) we present the mathematical model. In Section HI we 



analyze the stationary properties of the system and intro- 
duce a criterion based on the degree of non-orthogonality 
of the eigenmodes in order to identify the EP for finite 
system sizes. The beam dynamics is studied in Section 
|IV[ First, we study numerically the beam evolution in 
subsection |IV Aj while in the theoretical analysis is per- 
formed in subsection IV B Our conclusions are given at 
the last Section Iv] 



II. PHOTONIC LATTICE MODEL 

We consider a two dimensional honeycomb photonic 
lattice of coupled optical waveguides. Each waveguide 
supports only one mode, while light is transferred from 
waveguide to waveguide through optical tunneling. A 
schematic of the set-up is shown in Fig. [T] The lattice 
consist of two types of waveguides: type (A) made from 
lossy material (green) whereas type (B) exhibits the equal 
amount of gain (red) . Their arrangement in space is such 
that they form coupled (A-B) dimers with inter and intra- 
dimer couplings ta and t respectively. Such structure, 
apart from a global VT symmetry, respects also another 
anti-linear symmetry (in Ref. [33^ we coined this VdT- 
symmetry) which is related with the local T'T-symmetry 
of each individual dimer. 

Without loss of generality we rescale everything in 
units of the inter-dimer coupling ta = 1. In the 
tight binding description, the diffraction dynamics of the 
mode electric field amplitude ^'n.m = (in.m, ^n.m)"^ at 
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the (n, m)— th dimer evolves according to the following 
Schrodinger-like equation 



ear, 



K-i, 



■ db„ 



m+1 



+ tb„ 







dz 



H- ^'^bn,rn + fin+l,m+l + ^n+l,m-l + tan^m — 

(1) 

where e = + is related to the refractive index . 
Without loss of generality, we assume cq = and 7 > 0. 

For latter purposes, it is useful to work in the momen- 
tum space. To this end, we write the field amplitudes 
^„ m in their Fourier representation i.e. 

an,m{z) = J^^dka;dkyak^,ky{z)e-iip{i[nkx + mky]) 

bn,rniz) = J^^ dk^dkybk^ ^kyiz) e:xp{i[nkx + mky]) . 

Substitution of Eq. (l2| to Eq. (II]) leads to the following 
set of coupled differential equations in the momentum 
space 



where 



and 



d_ f d\^{z) 
'dz \ bi^{z) 



Hk = 



Dik) = -it + 2e 



ak{z) 
k{z) 



(3) 



-^7 D(k) 
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cos fcy); k : 



( kx , ky 



(5) 



In other words, because of the translational invariance of 
the system, the equations of motion in the Fourier repre- 
sentation break up into 2x2 blocks, one for each value 
of momentum k. The two-component wave functions for 
different k values are then decoupled, thus allowing for a 
simple theoretical description of the system. 



III. EIGENMODES ANALYSIS 

We start our analysis with the study of the stationary 
solutions corresponding to the Eq. ([s]). Substituting the 
stationary form 



idn,m, bn,mf = exp(-if B)^ 

in Eq. ^ we get 



(6) 



(7) 



The spectrum is obtained by requesting a non-trivial so- 
lution i.e. {A, B) 7^ 0. The corresponding dispersion 
relation [TUj has the form 



(8) 



For 7 = the dispersion relation is 

£ = ±\Dm (9) 



and we have two bands of width t 
pairs of diabolic points. 



2. There are three 



k^'^ = (±7r, ± arccos |) 
kj" = (0, ±(7r — arccos |)) 



(10) 



Expansion of -D(k) up to the first order around the DPs 
leads to 



(11) 



Substituting the above expression in the energy disper- 
sion given by Eqn. (|8]), we get the following linear rela- 
tion 



±{t^^i + {A-t^Wy) 



2N1/2 



(12) 



where "qx 



(fco. 



The standard passive (7 = 0) honeycomb lattice (zero 
strain) corresponds to t = 1. In this case, there are three 
pairs of DPs at (see Figj2]-a) 



k±'^ = (±7r,±f) 
k± = (0,±f ) 



(13) 



(4) For 1 < t < 2 the two pairs of DPs at kg start mov- 



ing toward each other while the pair at kj^ moves away 
from one another. At i = 2 a degeneracy occurs i.e. 
ki^'^ (±7r,0). At the same time the dispersion re- 



lation £ around kj'^ and is linear only in the kx 
direction (and quadratic in the ky). For t > 2 the two 
energy surfaces move away from each other and a gap 
between them is created (see Fig. [2]-b). Therefore for 
t > 2 the DPs disappeared for 7 = 0, and the conical 
diffraction is destroyed [M] ■ 

By introducing gain/loss to the system described by 
Eq. ([1]) the resulting effective Hamiltonian which de- 
scribes the paraxial evolution becomes non-Hermitian. 
In fact, for 1 < t < 2, any value of 7 results in 
complex eigenvalues, i.e. the system is in the broken 
T^T-symmetry phase. The resulting dispersion relation, 
resembles the dispersion relation of relativistic parti- 
cles with imaginary mass as was recently discussed in 
Ref. [ID]. 

In the case of t > 2 the size of the gap between the two 
bands can be controlled by manipulating the gain/loss 
parameter 7. In this case, there is a 7 domain, corre- 
sponding to the exact phase, for which the energies are 
real. It turns out from Eq. ([8| that the line 



IVT = t-2 



(14) 



defines the phase transition from exact to broken VT- 
symmetry jlOj . The mechanism for this symmetry break- 
ing is the crossing between levels, associated to the ex- 
ceptional points 



k± 



EP 



(0,±7r) 



k±p = (±7r,0) 



(15) 
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FIG. 2: Energy surfaces for a honeycomb lattice with ta ~ 1 
and (a) (t, 7) — (1,0) where a DP is present; (b) (t,7) = 
(2.5,0) where a band-gap has been created and (c) (t, 7) = 
(2.5, 0.5) where a gain/loss parameter led to the creation of a 
EP. 



(see Fig[2]-c) and belonging to different bands [33]: it 
follows from Eq. Q, that when 7 = j-pTi the gap disap- 
pears and the two (real) levels at the " inner" band-edges 
become degenerate. Evaluation of D{k^,ky) to second 



order in {r]x,riy) around the degeneracy points, leads to 

D{Tj,,7jy)^~i^vr + 2iv. + v^); v^ = vl + vl (16) 

thus resulting in the following dispersion relation 

(17) 



£ = ±./27prr?2 + (27pr + 4)7^2 



For large 7737-- values (i.e. 7737- ^ 2) one can approximate 
the above equation to get 



£ « ±v/27^77. 



(18) 



This comment will become important in the analysis 
of optical beam propagation discussed in the next section 

El 

Next, we turn to the analysis and characterization 
of the bi-orthogonal set of eigenvectors of our non- 
Hcrmitian system. The target here is to identify the prox- 
imity to the exceptional point in the case of finite Hilbert 
spaces, where finite size effects might play an important 
role in the analysis of the dynamics. The latter do not re- 
spect the standard (Euclidian) ortho-normalization con- 
dition. Let {Ln \ and |i?„) denote the left and right eigen- 
vectors of the non-hermitian Hamiltonian H correspond- 
ing to the eigenvalue i.e. 



{Lnl'H — {Ln\£n 

and 

^l-Rn) — £n\Rn}- 

The vectors can be normalized to satisfy 

(-^n|-^m) — ^nm 

while 



(19) 



(20) 



(21) 



Above is the dimension of the Hilbert space. 

An observable that measures the non-orthogonality of 
the modes, and can be used in order to identify the prox- 
imity to the EP in the presence of finite size effects, is 
the so-called Petermann factor which is defined as |35) 
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FIG. 3: Inverse Petermann factor for various system sizes (a) 
JV = 24; (b) JV = 168; (c) A/" = 440 . The asymptotic hne 
7737- — t — 2 (black hne) is approached by our numerical data 
(zero of K ) for increasing Af. The green line is plotted in 
order to guide the eye. The insets show representative K 's 
vs. 7 for a fixed coupling t = 2.5. 



We have studied the mean (diagonal) Petermann factor 



(23) 



n=l 



which takes the value 1 if the eigenfunctions of the sys- 
tem are orthogonal while is larger than one in the op- 
posite case. At the EP, a pair of eigenvectors associated 
with the corresponding degenerate eigenvalues coalesce, 
leading to a collapse of the Hilbert space. At this point, 
the Petermann factor diverges as 



l/l7-7-prl 



(24) 



K„ 



{Ln\Lm) {Rm\Rn) 



(22) 



[20I l36] . This indicates strong correlations between the 
spectrum and the eigenvectors which can affect drasti- 
cally the dynamics as we will see later. 

In Figs. [3^-c we present our calculations for K for 
different (t, 7) values and various system sizes. We find 
that as M increases the divergence is approaching the line 
77,7- — t~2, which was derived previously (see Eq. ( 14) 
above) for the case of infinite graphene. 



IV. DYNAMICS 

Armed with the previous knowledge about the eigen- 
modes properties of the T^T-symmetric graphene, we are 
now ready to study beam propagation in T^T-symmetric 
honeycomb lattices at the vicinity of the EPs. The ques- 
tion at hand is whether the collapse of the Hilbert space 
at the EP, will affect the CD pattern, and if yes, what is 
the emerging dynamical picture. 



A. Numerical Analysis 

We first study wave propagation in the honeycomb lat- 
tice numerically (Fig. |4]), by launching a beam with the 
structure of a Bloch mode associated with the EP, mul- 
tiplied by a Gaussian envelope. The Bloch modes at the 
tip can be constructed from pairs of plane waves with 
k vectors of opposite pairs of exceptional points. Thus, 
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interfering two plane- waves at angles associated with op- 
posite EP yields the phase structure of the modes from 
these points. Multiplying these waves by an envelope 
yields a superposition of Bloch modes in a region around 
these points. Figure |4] shows an example of the propaga- 
tion of a beam constructed to excite a Gaussian superpo- 
sition of Bloch modes around an EP. The input beam has 
a bell-shape structure, which, after some distance, trans- 
forms into the ring-like characteristic of conical diffrac- 
tion |34j . From there on, the ring is propagating in the 
lattice by keeping its width constant while its radius is 
increasing linearly with distance. The invariance of the 
ring thickness and structure manifests a (quasi-) linear 
dispersion relation above and below the EP (see Fig. 2); 
hence, the diffraction coefficient for wave packets con- 
structed from Bloch modes in that region is zero (infinite 
effective mass) . This is especially interesting because the 
ring itself is a manifestation of the dispersion properties 
at the EP itself, where the diffraction coefficient is in- 
finite (zero-effective mass). As a result, the ring forms 
a light cone in the lattice. The appearance of CD in 
the case of PT-lattices where the eigenvectors are non- 
orthogonal and coalesce at the EP singularity, provides a 
clear indication that the phenomenon is insensitive to the 
eigenmodes structure and it depends only on the proper- 
ties of the dispersion relation. 

The PT-symmetric conical diffraction shows some 
unique characteristics with respect to the CD found in 
the case of beam propagation around DPs for passive 
honeycomb lattices (i.e. 7 = 0). A profound difference is 
associated with the fact that now, the transverse speed 
of the cone is increased jTD] and in fact it can be con- 
trolled by the magnitude of the gain/loss parameter at 
the symmetry breaking point ^vT- This is shown in Fig. 
|4] where we compare the spreading of a CD for two dif- 
ferent 7737--values. 

In Figure [5] we report in a double logarithmic plot, our 
numerical measurements for the transverse velocity of the 
spreading ring for various 7737- values. The best linear 
fitting to the numerical data shows that transverse speed 
of the cone is proportional to ^7^7-. This behavior, can 
be understood qualitatively by realizing that the group 
velocity near the EP's is 
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FIG. 4: Propagation of a Gaussian superposition of Bloch 
modes associated with the vicinity of an EP at k = (0, ±7r), 
in a VT -symmetric honeycomb lattice. Shown is the beam 
intensity at normalized propagation distances of (a,d) z — 0; 
(b,e)2: — 10; (c,f) z — 15, for "fvT = 1 and t = 3 (upper 
panels) and y-pj- — 2 and t — 4 (lower panels). The input 
bell-shaped beam transforms into a ring-like structure of light 
of a non- varying thickness. 
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FIG. 5: Transverse velocity of the spreading ring versus the 
j-PT- Numerical simulation approve theoretical prediction for 
the transverse speed of the CD, which is Vg oc ^fypr- 



around the EP while it decays exponentially away from 
it. Specifically we assume 



(25) 



(ak(0),&u(0))^ = e-^V^(l,0)^ 



(26) 



(see Eq. (171 and Eq.(18l). At the same time we find Next, we calculate the evolution matrix U 



that the resulting VT-coue is brighter with respect to the 
one found in passive lattices [S] (i.e. the field intensity of 
the conical wavefront is larger) . 



where Hk is given by Eq. (4|. After a straightforward 
algebra and using the fact that 



Hk^ X 1, 



(27) 



B. Theoretical Considerations 

It is possible to gain valuable insight into the features 
of ■PT-conical diffraction by considering the field evolu- 
tion in the momentum space. We consider for simplicity, 
an initial distribution (ak(0), &k(0))^ that is symmetric 



where 1 is the unity matrix, we get 

U = cos(z|£:|)l - i{sm{z\£\)/\£\)Uu (28) 



Equation (28) is the starting point of our analysis. Sub- 
stituting Eqs. ( |4|16[ \l7\ , we find that the evolving am- 
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plitude of the field (a„ &„_m) is 



1 = 1,: 



(-l)'i [z - JvtHiT', m, z, g)] + g 



(29) 



0(n, m, z, 5)^/^ 
(j>{n, m, z, 5)^/^ 



i = l,2 



where 

(P{n,m,z,g) = [g + {-lYizf+n^ /{2j-pT+i)+m^ /{2jvt)- 

(30) 

Although our simplified calculations are not able to cap- 
ture all the features of the propagating cone, the above 
expression encompasses the main characteristics of the 
conical diffraction that we have observed in our numerical 
simulations. At z = 0, Eq. (29 1 resembles a Lorentzian, 



which slowly transforms into a ring of light, whose radius 
expands linearly with z with velocity ^27^7-, while its 
thickness remains unchanged. At the same time the field 
intensity on the ring in the case of T'T-symmctric lattices 
is brighter than the one corresponding to passive honey- 
comb lattices (i.e. vs behavior respectively). 



V. SUMMARY AND CONCLUDING REMARKS 

We studied numerically and analytically the prop- 
agation of waves in T'T-honeycomb photonic lattices, 
demonstrating the existence of conical diffraction arising 



solely from the presence of a spontaneous 'PT'-symmetry 
breaking phase transition point. In spite the fact that 
the eigenvectors are non-orthogonal and there is a col- 
lapse of the Hilbert space at the EP, the emerging cone, 
is brighter and moves faster than the corresponding one 
of the passive structure. Although, the realization of such 
photonic structures is currently a challenging task, active 
electronic circuits, like the one proposed in Ref. [TH], can 
be proven useful alternatives that might allow us to in- 
vestigate experimentally, wave propagation in extended 
T'T-symmetric lattices. 

These findings raise several intriguing questions. For 
example, how does nonlinearity affect PT-symmetric 
conical diffraction? What is the effect of disorder [37]? 
Is this behavior generic for any system at the sponta- 
neously T'T-symmetry breaking point? These intriguing 
questions are universal, and relate to any field in which 
waves can propagate in a periodic potential. It is ex- 
pected that in active metamaterials such phenomena will 
be present and they may have specific technological im- 
portance. 
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